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1. Introduction 

Unitary and special unitary groups play an important role in various fields of physics. 
Several problems that arise in the context of these groups require to express them in 
terms of a set of real parameters. In general, such parameterizations are not unique 
in the sense that the considered groups can be parameterized in many different ways. 
With regard to the diversity of problems it is reasonable to have a repertoire of different 
parameterizations available, in order to be able to choose the most convenient one for a 
given problem. In Ref. [T] we recently introduced the composite parameterization of the 
unitary group U{d) as an alternative to the canonical parameterization U = exp{iH) 
via hermitian matrices H and those presented in Refs. [21 |3l H]. It was shown that 
our parameterization enables a simple identification of redundant parameters when 
it is applied to describing orthonormal bases, density matrices of arbitrary rank and 
subspaces. For all these objects we found representations containing the minimal 
number of parameters needed. Due to its concise notation, simple implementation 
and computational benefits it has already found widespread applications in research 
on lattice correlation functions |5], quantum nonlocality |6], genuine multipartite 
entanglement [3 El E] and quantum secret sharing [10]. 

The aim of the present paper is twofold. First, we adopt our concept that was used 
in Ref. [1] to obtain a novel parameterization of the special unitary group SU{d). The 
need of additional representations of this group is not only apparent because of its vital 
relevance in all kinds of fields involving quantum physics (see Ref. [11] for an overview), 
but moreover is given because the number of available parameterizations is relatively 
low compared to lA{d). Here, our parameterization is proposed as an alternative to the 
canonical parameterization U = exp{iH) via traceless hermitian matrices H and the 
generalized Euler angle parameterization [12] introduced by Tilma and Sudarshan. 

Second, we rigorously derive the normalized Haar measure in terms of the 
introduced parameters for both the unitary U (d) and the special unitary group SU {d) 
of arbitrary dimension. In form of an infinitesimal volume element of a group, the 
Haar measure contains all information about the distribution density in its parameter 
representation. This in return is essential for the capability of generating uniformly 
distributed random unitaries, density matrices and subspaces, as they are important in, 
for instance, Monte Carlo simulations [13j or quantum data hiding [TH (TH]. Explicit 
expressions of the normalized Haar measure can furthermore be useful to tackle group 
integrals as they appear in lattice QCD [16], quantum optics [17], stochastic processes 
|18] and mesoscopic systems [T9]. In quantum information they can be found in the 
context of symmetric states {Werner and isotropic states [23 Ell [22]) and the a priori 
entanglement of quantum systems [221 1211 12S]- 

Within recent years much attention has been paid to integrals over unitary groups, 
i.e. /^(^) fiU, U*)dU whose integrand is a polynomial in U and U*. It was shown that 
such integrals can be replaced by a sum of function values f{Ui, U*) using a finite set of 
unitaries {Ui}^- Such sets are termed unitary t-designs [231261122], wherein t denotes 



Composite parameterization and Haar measure for all unitary and special unitary groups?) 



that the degree of the polynomial in U and U* is at most t. Whereas the existence 
of a unitary t-design was proven for all dimensions d and all polynomial degrees t (see 
Ref. |28j), it is generally unknown how to construct them for arbitrary d and t. In 
addition, there currently exists no simple analytic method for solving any integral over 
any polynomial. Due to this, there have been several attempts to find schemes to 
approximate unitary t-designs. Here, the usefulness of our results in the context of 
unitary t-designs and their approximations is obvious: Since we do not only provide the 
Haar measure but also the exact parameter ranges for group covering, the integration 
of polynomials can be performed explicitly and solved analytically in many cases. This 
in return allows to verify (or falsify) suggested unitary designs and to test the accuracy 
of approximations. Apart from that, our tools for computing integrals are not limited 
to polynomials, but are applicable to arbitrary functions. 

The paper is organized as follows. In Section [2] we review the composite 
parameterization of the unitary group U{d). In Section E] we introduce the composite 
parameterization of the special unitary group SU{d). In Section H] the general formula 
for the Haar measure on U{d) is stated and proven. An analogous formula for the Haar 
measure on the special unitary group SU (d) can be found in Section [51 Finally, in 
Section [6] we make useful remarks on computing integrals over U{d) and SU{d) using 
the composite parameterization and the associated Haar measure. 



2. Composite parameterization of the unitary group U{d) 



Consider a d-dimensional {d > 2) complex Hilbert space H = spanned by the 
orthonormal basis , . . . , \d)} . On this space define d one-dimensional projectors 



Pi = 10 (/| 

and d{d — l)/2 anti-symmetric matrices 

Ym,n = -i \m) {n\ + i \n) (m| 



KKd 



1 < m < n < d 



(1) 



(2) 



each acting on a two-dimensional subspace spanned by \m) and \n). In our previous 
paper |T], using these operators we have proven the following theorem: 

Theorem 1. Any operator of the unitary group U{d) can be written as[§| 



c 



d-l 



Y\_exp{iPiXi^i) 



.1=1 



(3) 



Y\_ I n exp (iP„A„,m) exp (iF„,„Am 

_m=l \n=m+l 

using d"^ real parameters {Xm,n}m,n=i,...,d in the ranges Xm,n £ [0, 2n] for m > n and 
>^m,n e [O, |] form<n. 

The idea behind this construction was to compose the unitary group out of 
'elementary operations' such as rotations and phase shifts. Here, these operations are 

I These can be considered as generalizations of the Pauh matrix a-y 



The order of the product is 0^=1 Ai = Ai ■ A2 ■ ■ ■ A 



N 
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realized through the matrix exponential functions exp {lYm^no) (generates a rotation) 
and exp (iPia) (shifts a phase). To make an ansatz for an arbitrary unitary operator 
it is useful to think of unitaries as orthonormal basis transformations. In this way, the 
form on]) can be interpreted as one option of incorporating d global phase operations 

exp{iPiXu) (4) 

as well as d{d — l)/2 rotations followed by relative phase shifts 

J^m,n — GXp {iPn^n.m) GXp {lYrn.n'^rn.n) (5) 

that result from partitioning the Hilbert space into two-dimensional subspaces according 
to ([2]). That this d^ parameter set of unitary operators indeed covers the whole unitary 
group U{d) was proven in Ref. This was done by showing that for any U E U{d) 
there exists a Uc such that UqU = 1. 



3. Composite parameterization of the special unitary group SU{d) 

Unitary operators obeying det U = 1 constitute a subgroup of U (d) called the special 
unitary group SU{d). The first new result of this paper is that this group can be 
parameterized similarly as U{d) using the concept of the composite parameterization. 
First, let us note that as special unitary operators satisfy an additional constraint the 
special unitary group SIA (d) can be described by — 1 real parameters. As redundant 
parameters are undesirable we have to find a parameterization that contains exactly 
this number of parameters. Second, it is known that U = exp{iHa) is special unitary 
for all a G M only if H is hermitian and traceless. In the composite parameterizaton of 
the unitary group U (d) we have used matrix exponentials of Pi and Ym,n as defined in 
Eq. dl]) and Eq. The latter is already traceless but any one- dimensional projector 
Pi has Tr(P/) = 1. As the projectors Pi were used to create phase shifts it is clear that 
they have to be replaced by a set of diagonal traceless operators. For that we introduce 
the following operators 

Zm,n = 1^) (^1 ^ 1^) ('^l 1 < m < n < d . (6) 

These are possible generalizations of the diagonal Pauli matrix acting on the subspace 
spanned by |m) and |n). In (|5]) the operation exp (iP„A„,m) was used to generate a 
relative phase shift between the vector components \m) and \n). However, the same 
effect can also be achieved with exp {iZm,n^n,m) meaning that (E]) can be replaced by 

A}ri,n — exp {iZjYi,nXn,m) 6Xp (^iYjYi^nXm,n) • (7) 

It now remains to turn our attention to the last part of (|3]), i.e. nf=i Gxp(iPiA/,/) . Here, 
each exp{iPiXi^i) is regarded as a global phase operation on |/). Special unitarity implies 
that there are only d — 1 independent global phase operations instead of d for U{d). 
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There is no unique way how these operations can be reahzed using matrix exponentials 
of ([n]). However, a possible and convenient choice is 



exp(iZi,rfA, 



Ll 



KKd-l . 



(8) 



In this version, the first d — 1 vectors |1) , . . . , |(i — 1) experience the phase shifts 



AX,, 



\d — 1) while the last vector \d) gets phase shifted in the overall 



inverse direction, i.e. e^*^i=i \d). Note that according to our labeling there is no 
parameter \d,d- Thus, in total we have the desired number oi d^ — 1 parameters \m,n- 
In summary, this leads to the next theorem: 

Theorem 2. Any operator of the special unitary group SU{d) can he written as 



c 



d-i 



n n exp (zZ„,„A„,„) exp (iFm,„A„,n) 



,m=l \n=m+l 



i-1 



]^exp(iZ/,rfA;,i) 



.1=1 



(9) 



using rf^ — 1 real parameters {Xm,n^ in the ranges Xm,n G [0, vr] for m > n, Am,n £ [O, |] 
for m < n and Xm,n ^ [0, 27r] for m = n. 

Proof. Proving that any U G SU{d) may be written as ([9]) is equivalent to showing 
that [/J,f/ = 1 can be achieved for all group elements. Let U = X]rs=i ^r,s k) {s\ be an 
arbitrary special unitary operator, i.e. U fulfils Ylt=i Kn fin,i = Yl'i=i '^im^i,n = ^mn and 
det U = 1. The conjugate transpose of Uc as given in ^ using the abbreviation ([7]) is 



U, 



c 



d-l 

JJexp( 



-iZd-i^dXd-i,d-i) 



1=1 



d-l 



n Ha 

m=l \n=l 



t 

d—m,d+l—n 



(10) 



The order of the factors in Uq implies that A\ 2 acts first on U. For U' = A\2U 
r) {s\ one obtains 



Sr,s=l ^r,s 



I 



e-*^2'^ cos(Ai,2)ai,, - e'^'''^ sin(Ai,2)a2 



i\2: 



''2,8 



e-*^2'i sin(Ai,2)ai,s + e*^^-^ cos(Ai,2)a2,s • 



i\2-_ 



(11) 

(12) 



a^^s for r > 2 if d > 2. We observe 



All other components remain unchanged, i.e. a'^^ 
that always be made zero using particular values for Ai_2 and A2.1: In case 

ai,i and a2,i both are zero, both parameters Ai,2 and A2,i can be chosen freely. If only 
ai^i = one chooses Ai_2 = |. If both are unequal zero then i vanishes for 



arg(e*^^2,iQ2 = arg(-ai,i) 

|02,l| 



tan(Ai_2 



(13) 
(14) 



This can always be achieved with A2,i G [0,7r] and Ai,2 G [0,f]- Analogously, one can 
make the component a'g ^ of U" = A| ^U' zero for the case d > 2. In this way, all 



Note that the indices m and n again run from 1 to d except that there is no Xd,d- 
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components ar^s with r > s of f/ = Yl 



d-l 
m=l 



T-\m At 

1 in=l ^ ^d-m,d+l-n 



u 



E 



d - 

r,s=l ^r,s 



\r) 



can be made zero. The unitarity constraints X]i=i ^ 'm i^n,i 



Ed —if — r 
i=l ^i,m^i,n ~ '-'mn 



imply that also all ar^s with r < s have become zero during this procedure. Hence, 
U is diagonal U = J2r=i'^r,r \r) {r\ where complex numbers of magnitude 

one, i.e. ar,r = e*"''. The first d — 1 phases ai, . . . ,ad-i can be compensated 

: ar- This is guaranteed to be 



via 



Ylf=l exp{-iZd-i4Xd-i,d-i) U by choosing A^, 



achievable with Ar.,r G [0, 27r]. Recall that so far we have only multiplied special unitary 
operators implying that UljU 



U is still a member of 



]\i=i exp{-iZd-i,dXd-i,d-i) 
SU{d). We have achieved that U^U is diagonal and that the first d — 1 diagonal entries 
are all equal 1. Now, the fact that det{UQU) = 1 still holds implies that also the last 
diagonal entry equals 1. Hence, Ul,U = 1, which proves the theorem. □ 



3.1. Remarks on the composite parameterzation 

It is worth mentioning some properties of the composite parameterization. In our 
previous paper pQ on the parameterization of U{d) we have shown that it can be rather 
insightful to gather the parameters Xrn,n in a matrix 



relative phases — )■ 



A, 



d,l 



A 



d,d 



^ rotations 



(15) 



In this notation the lower left entries represent relative phase shifts, the diagonal global 
phase shifts and the upper right rotations (with respect to the computational basis 
, . . . , \d)}) . Using this representation is particularly useful for illustrating which 
parameters are irrelevant for certain tasks. For instance, we have shown that for 
parameterizing an orthonormal set of k vectors {l^&i) , . . . , l^fc)} only the k{2d — k — 1) 
non-zero parameters of the following matrix are relevant 



Ai,2 

•^2,1 

Xk+i,i 
Xd,i 



A, 



d,k 



Ai 



fc+i 



Afc^fc+i 
Afc+i,fc 



Ai,d 

Xk,d 




(16) 



d ~ k 



d-k 



A similar example was given for parameterizing the set of fc-dimensional subspaces in 
where the corresponding 2k{d — k) relevant parameters Xm,n were illustrated. For the 
composite parameterization of the special unitary group SU{d) all this remains valid as 
the operations ([5]) and ([7]) are equivalent up to a global phase and are applied in the 
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same order in both cases. Using this matrix representation for the special unitary group 
one should only keep in mind that there is no diagonal element A^^^. 

Let us illustrate another important feature of the composite parameterization. The 
{d — kY oi {d — k)"^ — 1, respectively, non-zero parameters 





■ 

• 







Afc+i,fc+i 

^d,k+l 






Afc+i,(i 



A, 



d,d 



[17) 



d — k 



d-k 



correspond to the (special) unitary group for the {d — fc)-dimensional subspace defined 
by span(|fc + 1) , . . . , |(i)). This directly follows from the correct number of required 
parameters in combination with the fact that the subspace s]iaxi{Uc |1) , . . . , f/c \k)) is 
independent of the illustrated parameters. Note that this feature will be helpful in an 
upcoming proof. 



4. Haar measure on the unitary group U{d) 

Let us now assign an infinitesimal volume element dUd to the unitary group U{d) in 
terms of the composite parameterization Uc = Uc{\i,i, ■ ■ ■ , Xd,d)- This can be achieved 
by determining the associated Haar measure. This means that, as for any compact 
Lie group, we must find a measure of volume which is left and right invariant [291 130] . 
Explicitly, we require that dUd satisfies 



dUd = d{Uc) = d{U,Uc) = d{UcU2) 



2) , 



for all Ui,U2 € U{d). Generally, an invariant measure is (up to an irrelevant constant) 
determined by the absolute value of the Jacobian determinant 



Jd 



|det(jfcj)| 



det 



d{ai, 



wherein {uk} are coefficients of the unitary U 
orthogonal operator basis {bk} of x C^, i.e. 



U{ai 



Uk 



Mbibk) 



U 




(19) 

arf2 ) expanded in an 
(20) 



and where {a;} are any d'^ parameters that cover U{d). Here, a left or right translation 
Ui, U2 G U{d) merely induces a unitary basis transformation of {bk} which is length and 
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angle preserving. Hence, Jd is invariant under these transformations and 

dUd = Jdl[dai (21) 

1=1 

is a Haar measur^. Our aim is to derive a general expression of dUd for arbitrary d in 
terms of the parameterization introduced in Section [21 i.e. 

d 

dUd = ^Y\ d\k,i , (22) 

where Nd is a normalization constant such that /^^^^ dUd = 1- We obtain the following: 

Theorem 3. In terms of the d"^ parameters Xm,n introduced in Theorem U\ the 
(normalized) Haar measure on the unitary group U{d) reads 

d—l d d 

n ™i(A„,„)cos"'"-"'''(A„,„) J] dAy . (23) 



m=l n=m+l 



with 



- (^')^'^"'"' (24) 

llm=llln=m+l2(^-"^) 

such that /^(^-) dUd = 1. 

Proof. For simplicity, let us start with the case d = 2. Using the canonical operator 
basis {bk} = {|1) (1| , |1) (2| , |2) (1| , |2) (2|} and the order of the parameters = 
{^1,1, Ai,27 ^2,1, A2,2} we obtain the Jacobian matrix 

d{Xl,l, Ai,2, ^2,1, A2,2) 

/ ze^^i.i cosAi,2 -e^^i-i sinAi,2 \ 

e*^2,2 cos Ai,2 ie'^^'^sinXi 

y -e^^2.i+iA2.2 gin;^^ 2 ze^^2,i+»A2,2 cog_x^ 2 26*^2,1+^2,2 ^^g ^^^ ^ y 

Using the Laplace expansion and elementary simplifications one finds 

J2 = 2 sin(Ai_2) cos(Ai^2) = sin(2Ai^2) • (26) 
The relation ( 122|) combined with normalization 

"27r /.27r /.7r/2 /■27r 7 

/ / / J2 



/ 



/ / / Tf ^Ai, irfAi,2ciA2,irfA2,2 = 1 , (27) 

A2 2=0 JX2 1=0 ~'Ai.2=0 JA,.i=0 ^^2 



^ Intuitively: When changing from the orthogonal coordinates {ufc} to the non-orthogonal coordinates 

3(mi,...,Mj2) 

d(ai,....a^2) 



{ai} the infinitesimal volume element transforms as Y[k=i^^k — det g^q^' "'q''^) Ilf^i according 



to the Jacobian determinant. 
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yields the normalized Haar measure 

1 



47]-3 



sin(Ai,2) cos(Ai,2)rfAi,i(iAi_2(iA2,ic?A2,2 



f28) 



which is in accordance with Theorem [31 

One could in principle compute dUd analogously for arbitrary d, i.e. using 
the canonical operator basis {6^} = {|1) (1| , |1) (2| , |1) (3| , . . . , \d) (rf — 1| , \d) {d\} and 
the naive order {ai} = {Ai,i, Ai,2, Ai_3, . . . , Ad,d_i, A^^^i}. For instance, a long and 
cumbersome computation reveals that 

1 ^ 
dUs = — sin(Ai,2) cos(Ai,2) sin(Ai,3) cos^(Ai,3) sin(A2,3) cos(A2,3) Y\. ^-^^.^ (^9) 



kA=l 



demonstrating the validity of Theorem [3] for d = 3. Unfortunately, in this way, 
computing the determinant of the d'^ x d"^ Jacobian matrix becomes increasingly 
unfeasible the larger d gets. More importantly, this approach is not suitable to prove 
a general expression such as (!23|) for all d. However, the Jacobian matrix can be 
considerably simplified by taking into account the invariance of the Haar measure, the 
structure of the composite parameterization as well as the freedom in the choice of the 
operator basis and the order of the derivatives d/dX^^y In this way, the correctness of 
Theorem [3] can be verified for all d. 

First, due to the left invariance of the Jacobian determinant Jd = |det(jfc^;)| we are 
allowed to perform any transformation 



Jk,i 



duk _ 

dai ~ Tiiblbk) 



Jk.l 



(30) 



t dUc 
C dXx.v 



with Ui G U{d). Here, it is beneficial to choose Ui = —iUl; since the matrix —iU, 

has a simpler form due to the fact that Uq and ■§^^ cancel each other out partially. For 
instance, since all projectors Pi = \l) {l\ commute, for any derivative with respect to a 
global phase transformation d/dXi^i one obtains 



-tU. 



^9A 



Ll 



-lU^UciPi = Pi = \l) {l\ ■ 



(31) 



As can directly be inferred from the structure of Uc ([3]), the derivatives dUc/dX^^y are 



x-l 



n n 

_m=l n=m+l 

for X < y, and 

'y-l d 

n n A, 

m=l n=m+l 



n A. 



_n=x+l 



i^x,\ 



n A. 



n=y+l 



d-l 



n n A, 



■ x-l 

n ^y^^ 

n=y+l 



tPx 



J Lm=x+1 n=m+l 



■ d-l d 

n n ^ " 

m=y+l n=m+l 



]^exp(iPiA«,z) 



1=1 



Y\_exp{iPiXi^i) 



1=1 
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for X > y. Consequently, if we apply — «f/^ fro^^ the left, the products to the left of iY^^y 
(respectively iP^) cancel out and we get 



W P U 

^x,y-^ x'^x,y 



x,y 



for X < y , 
for X > y , 



(32) 



where 



Ux,y 



n=y+l ^2),n 



n=x ^y,n 



yrd-1 yrd 
1 Lm=x+1 1 Ir 



A. 



m=x+l 1 lra=r?i+l »Ti,n 

yrd 
1 Lm=y+1 i i 



n=m+l Am,™ 



n?=xexp(«P/A,,0 
n?=« exp(iPA/,/) 



for X < y 
for X > y 



(33) 



Here, it is important to realize that UlyYx^yUx,y and UlyPxUx,y do no longer contain 
operations Am,,n with m < min{x,?/}, meaning that there are no off-diagonal elements 
|m) (n| and \n) (m| with m < min{x, y}. This and the observation (13T]) imply that when 
we choose the following orthogonal operator basis and order 



hi 


|1)(2| + |2)(1| , 


«! 


= -^2,1 , 




-z|l)(2|+^|2)(l| , 


0:2 


= ^1,2 , 


h 


|1)(3| + |3)(1| , 


aa 


= ^3,1 5 




-z|l)(3|+^|3)(l| , 


0:4 


= ^1,3 ) 


&2(ci-l)-l = 


|1) {d\ + \d) (1| , 


a2(d-i)-i 


= ^d,l , 


&2(d-l) = 


{d\+i\d) {1\ , 


Oi2{d~l) 


= "^14 , 


^2(d-l)+l = 


|2) (3| + |3) (2| , 


«2(d-l)+l 


= -^3,2 , 


^2(ci-l)+2 = 


-z|2)(3|+«|3)(2| , 


«2(d-l)+2 


= -^2,3 5 


^2(d-l)+3 = 


|2) (4| + |4) (2| , 


«2(d-l)+3 


= -^4^2 ; 


^2((i-l)+4 = 


-^|2)(4|+^|4)(2| , 


tt2(cZ-l)+4 


= -^2,4 , 


bd^-d-i = 


\d - 1) {d\ + \d) {d - 1\ , 


ad'^-d-i 


= ^d,d-l 


bd-2-d = 


-i\d-l) {d\+i\d) {d-l\ , 


ad^-d 


= ^d-l,d 


bd^-d+i = 




Oid^-d+l 


= ^1,1 5 


bd'^-d+2 = 


|2) (2| , 


ttd2-d+2 


= -^2,2 J 


bd^ = 




arf2 


= ^d,d , 



(34) 



the Jacobian matrix becomes a lower block-triangular matrix 
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if, 



kl) 




M,.2 
4x4 








Md-, 
2x2 






D 
dxd 



where all entries outside the grey-shaded blocks are zero. Thus, the Jacobian 
determinant simplifies to a product of the determinants of the blocks Mj and D 

Jd = |det(^)| = [\[ IdetM.I j |detD| . (35) 

Furthermore, we have that |detD| = 1 since D is a c? x d identity matrix due to fl31l) 
and the choice (IMI) . It now remains to investigate the blocks Mj. Let us first consider 



-iU, 



t dUc 



for the Hilbert space "H' = C whose dimension is lower by one. In this 



case we have 



-iU, 



dU, 



c 



c 



dX 



x,y 



where 



x,y 



Y[n=y+1 



U'L=x+l UnJm+l ^m,n Utx exp(iPA 



W Y U 
iP P u 

^ x,y-^ x^x,y 



id-l 



for X < y 
for X > y 



(36) 



nn=x ^y,n Y[ni=v+1 11 



ln=m+l ^rn,n 



nt^exp (iPA/,z) 



for X < y 
for X > y . 



(37) 



By substituting A^.n — ^ Am.+i,n+i, \m) — )■ |m + 1) and (m| — )• (m + 1| one can directly 
see that these matrices are equal to —iUj^-ir^ with x,y > 2 and Uc for the Hilbert 
space "H = of full dimension d. Hence, since = |detMi| Y[iZ2 |detMi| we have 
established the recursion formula 

Jd = IdetMil Jd-i , 
where the parameters in J^-i are to be substituted according to Xm n ^ X 



(38) 

m+l,ra+l- This 

relation is a direct consequence of the fact that each —iUt-^r^ with x,y > 2 only 
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contains parameters Xm,n with m,n > 2. As discussed in Section l3TT] the set of {d — 1)^ 
parameters Xm,n satisfying m,n > 2 correspond to the unitary group U{d — 1) for the 
{d — l)-dimensional subspace span(|2) , . . . , \d)). Consequently, (138|) simply reflects the 
well-known fact [121 12] that the infinitesimal volume element dUd on U{d) is the product 
of the infinitesimal volume element dUd-i on U{d—1) times a to-be-determined function 
g{d), i.e. 

dUd = g{d)dUd-i . (39) 

Now, due to the above mentioned reasons, for the composite parameterization dUd-i is 
a function of the parameters \m,n with m,n>2. In addition to that, g(d) is a function 
of the 2d — 1 parameters {Ai^i, Ai^25 • • • , Xi,d} and {A2,i, . . . , A^^i} since they extend the 
unitary group — acting on span(|2) , . . . , \d)) to the entire U{d) on "H = C^, hence, 
g{d) = g{d; Ai^i, Ai_2, • • • , Ai^^; A2,i, . . . , A^^i). For the sake of clarity, let us illustrate this 
by the example dU2 — > dU^ by comparing (125]1 and dU2 is (up to a constant) given 
by sin Ai,2 cos Ai,2C?Ai,i(iAi,2C?A2,i(iA2,2- According to our considerations dU^ must contain 
the same expression in terms of A2,25 A2,3, A3^25 A3 3 and it must look like 

dUs ~ g{3; Ai,i, Ai,2, Ai,3; A2,i, A3,i) sin A2,3 cos A2,3'^A2,2C?A2,3C?A3,2ciA3,3 . (40) 

Equation fl29|) shows that this is indeed the case since it can easily be observed that 

g{3; Ai,i, Ai,2, Ai,3; A2,i, A3,i) ~ sin Ai,2 cos Ai,2 sin Ai,3 cos^ Ai,3dAi,i(iAi,2rfAi,3(iA2,irfA3,i . 

It now remains to find a general expression for |det Mi \ for arbitrary d, since the formula 
of Jd can then be derive via (155]) and induction. In general, the matrix Mi does not have 
a simple form which makes it difficult to compute its determinant. Here, it is useful to 
exploit that, first, there is no explicit dependence of Jd on parameters corresponding 
to phase operations, i.e. all Xm,n with m > n. This can directly be followed from the 
construction of the composite parameterization in combination with the Haar measure 
on U{2) and the Jacobian matrix fl2^ . Since in the Jacobian matrix fl25p and in 
its determinant the phase operations Ai,i, A2,i, A2,2 only appear as e*'^^'^ e*^^'\ e*^^'^, 
i.e. complex numbers of magnitude 1, and since in the end we are only interested 
in the absolute value of this determinant, it is clear that dU2 ( l28|) is independent of 
these parameters. As the composite parameterization only combines W(2)-operations 
by incorporating all possible 2-dimensional subspaces in the same statement holds of 
course true for all d (see for instance fl29|l for the case d = 3). Second, the circumstance 
that g{d) is independent of all parameters Xm,n with mm{m,n} > 2 implies that these 
parameters do not affect |detMi|E|. For these reasons one can set the parameters Xm,n 
with min{m, n} > 2 and m > n to zero in fl33p without altering | detMi|, i.e. 



detMil 



det (Mi|^,^^„=o|Kn)^{(i,2)....,(M)}}) I = I detMil . (41) 



The independence of |det Mi \ on Xm.n with minim, n} > 2 can also be confirmed by exploiting again 
the left and right invariance of the Haar measure. 
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Here, fl32l) reduces to 



{A,„,„=0|{m,n)^{(l,2),...,(l,d)}} 



Y\ exp(iYi,„Ai,„) 

,n=3/+l 



t 



n=y+l 



1,2/ 



{A„,„=0|{m,n)^{(l,2),...,(l,d)}} 



]^exp(iFi,„Ai,„) 



n=x 



P. 



where O^ i and Oi^y are orthogonal (real) matrices since they are a product of the 
operations exp(zYi^„Ai,„) which explicitly read 

cos(Ai,„) |1) (1| + sin(Ai,„) |1) (n| - sin(Ai,„) \n) (1| + cos(Am,„) \n) (n| + ^ \k) {k\ . 

(42) 



According to fl34|) . the 2{d — 1) x 2{d — 1) elements of Mi are now determined by 



Tr {hlOl^Yi,yOi,y) = Tr (Oi,ybiOlyYi,y 



and 



Tr (6lOj,P.O,,i) = Tr {o^AoliP. 
with x,y ^ {2, . . . , c?}, and 



(43) 



(44) 



= bl 



|1) ((A; + 3)/2| + |(A; + 3)/2) (1| = Xi,(,+3)/2 - odd , 
-^|l)((A; + 2)/2|+^|(fc + 2)/2)(l| = ri,(fc+2)/2 - even 



(45) 



with k G {1, . . . ,2{d — 1)}. Now consider the coefficients and depending on 
different Xi,^ and ^i,^ (2 < m < d) of (gS]): 



• Tr (Oi^yXi^m,0^yyL^y) for arbitrary m: 

Tr (Oi,,Xi,^0?;,ri,,) = . 

Note: Oi^yXi^mOiy is symmetric while Yi^y is antisymmetric. 

• Tr [Oi^yYi^rnOjyYi^y) for m < 

Tr {Oi^yYi^^OlyYi^y) = . 

Note: Oi^yYi^rnOi y is orthogonal to Yi^y ioi m <y since the product in 
Oi,y = nn=s/+i6xp(«>^i,nAi,n) starts with n = y + l. 



(46) 



(47) 
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• Tr (Ox,iyi,mOjiPj;) for arbitrary m: 

IV (0.,iYi,^OjiP,.) = {x\ 0.,iYi,„.Oji \x) , (48) 

= -i {x\ 0,,i |1) (m| Oji |x) + i {x\ 0,,i |m) (1| Oj^ \x) , (49) 

= . (50) 

Note: 0x^1 is an orthogonal (real) matrix. 

• Tr (0,,iXi,^Oj,iP,) for m < x: 

{0,,,X^,mOl^P,) = {x\ O^^^X.^mOl^ \x) , (51) 

= (x| |1) (m| Oji |x) + {x\ |m) (1| Oj^ \x) , (52) 

= . (53) 



Note: Ox,i = nn=x ^^P(^^i,n-^i,n) off-diagonal elements {x\ Ox,i \m) for m < x. 

These four observations imply for the operator basis and parameter order 

h = |1)(2| + |2)(1| , 

62 = -z|l)(2|+^|2)(l| , 

63 = |1)(3| + |3)(1| , 

64 = -z|l)(3|+^|3)(l| , 

b2(d-iyi = \l){d\ + \d){l\ , 
b2(d-i) = -t\l){d\+i\d){l\ , 

that Ml is a lower triangular matrix 

/ 



— {jk,l)k,l = l,...,2(d-l) — 

For the determinant we can now restrict on determining the diagonal entries of Mi 
which are given by 

. Tr {Oi,,Yi,^Ol^Yi,y) I Tr{YlJ for m = y: 
IV (Oi,^yi,^0?;^yi,^) = (1| Oi^ruYi^mOlJ'i^ru |1) + (m| Oi^ruYi^mOlJ'i,^ \m) , (55) 

=i (1| Ol,m^l,mO?^m 1"^) " ^ ("^1 ^i.^^^i.^O?^^ |l) , (56) 

= (1| Oi,^ |1) (m| O?;^ |m) - (1| |m) (1| |m) 

-(m|Oi,^|l)(m|0?;^|l) + (m|Oi,^|m)(l|0?;^|l) , (57) 
=2(l|Oi,^|l)(m|Oi,^|m) . (58) 
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Note: Oi^m = nn=m+i ^^P(^^i."'^i:") ^'^^^ have off-diagonal element (l|Oi,m|"^) 
since the product starts with n = m + 1. 
• Tr (0,,iXi,^Oj^,P,.) /Tr(X2„) for m = x: 

Tr (Om,iXi,mO^ iPm) = (m| Om,iXi^rnO'^ i \m) , (59) 

= {rn\Ora^\l) {m\0'^^\rn) + {rn\Orn,i\m) {l\0'^^\rn) , (60) 



2 (m| Om,i |1) {m\ Om,i \m) 



(61) 



We have thus found simple relations between the diagonal entries of Mi and the matrix 
elements {2 < m < d) 

d 

(l|Oi,^|l)= J] cosAi,,, (62) 

n=m+l 

(m| Oi .m |m) = 1 , 

("^1 |1) = -sinAi,m W cosAi,„ , 

n=m+l 

(m| 0^,1 \m) = cosAi,^ , 

which can easily be obtained via the definitions of Om,i and Oi^m together with 
Hence, since Tr(Xf ^) = Tr(Y]^^) = 2, we find that the diagonal entries of Mi are 

d 

Tr(0i,™yi,„0?;^Fi,^)/Tr(ri2j= JJ cosAi,„ , 

ri=m+l 

Tr (Om,iXi,„0^ iPm) /Tr(Xi ,„) = -sinAi,™ Y\_ cosAi,n • 

n=m 

If we multiply all these entries {m = 2, . . . ,d) we finally obtain 



(63) 

(64) 
(65) 



IdetMil 



k=2 \n=k+l 



JJ sin Ai,, JJ COS Ai,„ 



.1=2 



nsin(Ai,„)cos2("-i)-i(Ai,„) . 



n=2 

Using the recursion formula (138|) and induction one finds the final result 

d-l d 

Jd=n n sin(A^,0cos2("-™^-'(A^,„) • 

m=l n—m+l 

The corresponding normalizing constant X^ is given by the integral 

d 

'"W ki=i 



/ Jd\\ dXk,i , 
Ju(d) 



2n {k>l) ptt/2 {k<l) 
Afc,i=0 J\k,i=0 



d-l d 



n n 

m=l n=m+l 



COS 



2(n— m) — 1 



(A. 



Yl dX 



k,l 



(66) 
(67) 

(68) 
(69) 

(70) 

(71) 
(72) 



fcj=i 
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which can straightforwardly be computed 

d-l d , tt/2 

_nos2('^-'^)rA^„^ 

2{n — m 

m=l n=m+l 



AT, =(27r)'^('^+i)/2 ff n . ^ . cos^(— )(A^,„) 



nliinl=m+i2(^-H 



(73) 
(74) 

□ 



5. Haar measure on the special unitary group SU{d) 

Theorem 4. In terms of the d^ — 1 parameters Xm,n introduced in Theorem 2 the 
(normalized) Haar measure on the special unitary group SlA{d) reads 

d-l d 

dUd = j^X{ n sin(A^,„) cos2("— (A^,„) J] d\k,i , (75) 

m=l n=m+l k,l 

with 

od-l d(a!+l)/2-l 

Nd = 2 (76) 

Um=l nn=m+l2(«-H 

such that /j^^^) dUd = 1. 

Proof. The proof is similar to the previous one — only minor modifications have to be 
made. Given the special unitary group SU{d) in parameterized form ?7(ai, . . . , ad2_i), 
to construct the Haar measure one must determine the absolute value Jd = |det(jfc,i)| of 
the determinant of the Jacobian matrix 

_ d{ui,...,Ud2.r) 
d[ai, . . . ,arf2_i) 

wherein {uk\ are coefficients of the special unitary operator U{ai, . . . , ad'^-i) expanded 
in an orthogonal basis {6^} of traceless operatorqj of x C^, i.e. 



Uk 



Ti{h\bk) \ dai ^ dai 

In this terminology, the (normalized) Haar measure reads 

J d^-i 

dUd = ^l[ • C^S) 

* A basis for the vector space of operators x has (P elements. The constraint det J7 = 1 on 
special unitary operators implies that the trace of any derivative dU/dai with respect to an arbitrary 
parameter is always zero. Traceless operators form a (P — 1 dimensional subspace of x C*. As we 
have (P — 1 parameters {a;} it is required to express the derivatives dU /dai in a basis of this subspace 
to obtain cP — 1 linearly independent column vectors. 
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Our aim is to derive a general expression of dUd for arbitrary d in terms of the 
parameterization introduced in Sectional i.e. 



dU, 



k,i 



As in the previous proof we make use of the left invariance of the Haar measure on 
Sl({d), i.e. Jd = I det(jfc,i)| = | det(jX,)| where 



-Iriblh 



Jk,l 



Jkl 



Mbibk) 



If for the composite parameterization (Theorem [2]) Ui is chosen to be —iUl< one obtains 
analogously to fl3Tl) - fl33|l that 



-tU. 



dU, 



c 



c 



x,y 



W Y U 

^x,y^ x,y'-^x,y 



^x,y^y,xUx,', 



for X < y 
for X = y 
for X > y 



(82) 



where 



a 



x,y 



n 



n=y+l 



Y\.m=x+i n 



X-+1 1 l7i=m+l 



A. 



n?=iexp(iZi,rfA 



for X < y 
for X > y . 



(83) 



These operators are now to be expressed in an orthogonal basis {bk} of traceless 
operators. Since the first d"^ — d operators that were used in (IMIl already are mutually 
orthogonal and traceless, only the d diagonal operators \k) {k\ with 1 < k < d have to 
be replaced by ci — 1 diagonal operators with vanishing trace. A convenient choice are 
the d — 1 mutually orthogonal operator^ 



{d - k){d - k + 1) 



id-k)\k){k\+ J2 1^) 



n\ 



(84) 



n=k+l 



Note that these are the generahzed diagonal Gell-Mann matrices |31j in reversed order. 
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with l<k<d-l obeying Tr{LlLk) = Tr{Ll) = 2. For the order 



h 


|1)(2| + |2)(1| , 


a-i — 


^^2,1 , 


62 


-^|1)(2|+.|2)(1| , 


a.2 — 


^^1,2 , 


63 


|1)(3| + |3)(1| , 


as = 


^^3,1 ) 


h 


-i|l)(3|+i|3)(l| , 


— 


^^1,3 , 




\l){d\ + \d){l\ , 


(^2{d-l)-l = 


^d,l , 


b2{d-i) = 


-i\l){d\+i\d){l\ , 


Q;2(d-i) = 


^l,d 1 


b2(d-i)+i = 


|2) (3| + |3) (2| , 


Oi2{d-l)+l = 


, 


^'2((i-l)+2 = 


-2|2)(3|+^|3)(2| , 


Q;2((i-i)+2 = 


^2,3 , 




|2\ (41 + |4\ (21 




■^4,2 ) 


^2(d-l)+4 = 


-i|2)(4|+i|4)(2| , 


«2(d-l)+4 = 


-^2,4 , 


bd^-d-l = 


\d - 1) {d\ + \d) {d - 1\ , 




^d,d-l 


bd^-d = 


-i|d-l)(ci|+i|d)(rf-l| , 


ad2-d = 


^d-l,d 


bd'2-d+i — 


^1 , 


o^d^-d+i — 


^1,1 ) 


bd'^-d+2 = 


L2, 


(^d'^-d+2 = 


''^2,2 , 


bd?-i = 


Ld-i , 


Oid?-l = 


^d-l,d-l 



^5) 



the {d^ — 1) X {d^ — 1) Jacobian matrix is again lower block-triangular 



ifk,i) — 




M,.2 
4x4 








2x2 






D 

d-1 X d-1 
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Here, the block D is not diagonal but lower triangular since 





Tr(L)jZm,d) 



for k < m 
^ -{k — d — 1) for k = m 

for k > m . 



-k){d-k+l) ' 



^6) 



(d-k){d-k+l) 

Since the diagonal entries of D are merely real numbers we find again that 



d-l 



Jd= |det(j^_;)| =QjJ|detMi 

i=l 



where q is a constant that can be dropped since the Haar measure will be normalized at 
the end anyhow. As the composite parameterization of the unitary group and the special 
unitary group share the same structure the discussion between (138|) - ( HTj) essentially 
remains the same, meaning that (ignoring irrelevant constants) one finds again the 
recursion formula 



where Mi = Mi 



Jd= |detMi| Jd_i , (88) 
{A„,„=o|(m,n)^{(i,2),...,{i,d)}} ^^^^ A^.n A„+i,„+i in Jd-i- Here, the 



relevant operators of ([H2D have the form 



{A^,„=0|(m,n)^{(l,2),...,(l,d)}} 



n t 



J]^ exp(zyL,„Ai,^ 

n=y+l 



{Xm,n=0\{m,n)(^{(l,2),...,{l,d)}} 



JJexp(zFi,„A 



l,nJ 



J]^ exp(zFi,„Ai, 

n=y+l 



JJexp(zYi,„Ai,^ 



containing the orthogonal matrices Ox,i and Oi^y that also appeared in the previous 
proof. The 2{d — 1) x 2{d — 1) elements of Mi are determined by (compare with the 
order 



Tr [h\Ol^Y,,yO,,y) = Tr (o,,yblOlyYi,y 



and 



Tr {blOl,Z,,xOx,i) = Tr (OxAol^Z.^x) 
with X, ?/ G {2, . . . , d}, and 



(89) 



(90) 



bk = bi 



|l)((fc + 3)/2| + |(fc + 3)/2)(l| 



^i,(fc+3)/2 k - odd 



-«|l)((A; + 2)/2|+2|(A; + 2)/2)(l| = Fi,(a.+2)/2 k 



even 



(91) 
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with k G {1, . . . , 2((i— 1)}. Now consider the coefficients of f l89|) and f l90|) for the different 
operator basis elements occurring in ( MT\i : First, notice that Mi is again lower triangular 
as in the previous proof since (189|) and (H3l) are identical and 

• Tr (O^^iYi^mOj^Zi = for arbitrary m. 

Note: Oi^yYi^rnOiy is antisymmetric, while Zi^x is symmetric. 

• Tr (0^,iXi^„OjiZi^^) = for m < x: 

Tr (0,,iXi,^OjiZi,,) = (1| 0,,iXi,^Oji |1) - (x| 0,,,X,,^Ol, \x) , (92) 

= (l|0.,i|l)HOjJl)-(l|0.,i|m) (IIOJJI) 

- (x| 0^,1 11) (m| Oj,i \x) - {x\ 0,,i |m) (1| O^, \x) , (93) 

=0 . (94) 

Note: Ox,i = nn=a; ^^P(^^i.n-^i,n) ^^oes uot havc off-diagonal elements (1| Oxa \m) and 
{x \ Ox A \^) for m < X. 

Thus, it again suffices to compute the diagonal entries of Mi, half of which are already 
known from fl66ll 



Ti{Oi,mYi,mOlJ'i,J) /Ti{YlJ= W cosAi,„. (95) 

n=m+l 

The remaining ones are found to be 

d 

Tr (Oi,^Xi,^0^„Zi,^) / Tr(Xj„) = 2 sin Ai,^ JJ cos Ai,„ , (96) 

n—m 

which directly follow from 

Tr (Oi,^Xi,„O^^Zi,„) = (1| 0„,iXi,^0^_i |1) - (m| 0„,iXi,^0^_i \m) , (97) 

= (1| 0„,i |1) (m| 0^,1 |1) + (1| Om,i \m) (1| 0^,1 11) 

- (m| Om,i |1) ("^1 05^,i 1"^) - ("^1 0mA \m) (1| O^^i \m) , 
=2 (1| Om,i \m) (1| Om,i |1) - 2 (m| 0^,1 |1) (m| Om,i \m) . 

(98) 

and the definition of 0^,1 together with fH2]) 

(1| Orri,i |m) = sin Ai,™ , (99) 
d 

(l|0™,i|l) = J] cosAi,, , (100) 

n=m 

d 

(m| Om,i |1) = -sinAi,m W cosAi,„, (101) 

n=m+l 

(m| Om,i \m) = cosAi,m ■ (102) 
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Since fl96|) differs only by the factor two from f l66|) we again obtain for Ja the result (jTO 

d-l d 



Jd = c'd JJ JJ sin(A„,„) cos' 



2{n-m)-l 



(A. 



(103) 



m=l n=m+l 

up to an irrelevant multiplicative constant c'^. Hence, the Haar measure on SU{d) reads 

d-l d 



^ d-l d 

dUd = ^l[ n sin(A^,„)cos2(-— )-i(A„,„)nrfA,,z 



(104) 



m=l n=m+l 

The normalizing constant is given by the integral 



2-K (k=l) j-TT (fc>0 /.7r/2 {k<l) 



Ati=0 JXk,=0 JXk.i=0 



d-l d 



n n 

,m=l n=m+l 



COS 



2(n-m)-l 



( Arra,7i) 



kl 



Since there are c? — 1 parameters A^,; with k = I and d{d — l)/2 parameters A^^^ with 
k > I one finally obtains 



-1 



■ 2(?7, — m) 

1 n=m+l ^ ' 
2d-lj^d(d-^l)l2-l 



COS 



2(n— m) 



( Am,n) 



7r/2 



(105) 

(106) 
□ 



6. Remarks on integrals over unitary groups 

As previously mentioned, our results can be used to compute group integrals, i.e. 
integrals of the form J f{U,U*)dU where one integrates over the entire group U{d) 
or SU{d), respectively. At least three things are needed when one intends to 
explicitly compute such integrals: A parameterization of the corresponding group, exact 
knowledge of the parameter ranges and the normalized Haar measure. All this is 
provided in the present paper. Theorems [T]-|l] can straightforwardly be applied without 
knowledge of further technicalities (e.g. details appearing in the proofs). Whether or 
not a given integral can be solved analytically in this way of course depends on the 
integrand. However, for many physical problems the function f{U, U*) is a polynomial in 
the components of U and U* . In this case, when U and dU are inserted in parameterized 
form according to Theorems [T]-|U we have that the integrand is a polynomial in cos Xm,n, 
sinAm.n and e^*'^"-". Neglecting the computational effort, such integrals can always 
be solved analytically (see Ref. [32] and references therein). Besides that, our results 
constitute a good starting point for the integration of non-polynomial functions using 
numerical methods. 

A detailed analysis on integrals that can be solved in this way shall be presented in a 
subsequent paper. We are convinced that due to the simplicity of the parameterization 
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and the associated Haar measure, it is possible to find several general results and to 
gain a better understanding of integrals over U{d) / SlA{d). In order not to go beyond 
the scope of this paper we conclude with simple examples that can be compared with 
existing results consistency check. 

6.1. Example 1 

In Ref. [33] it was shown that 



I (l|f/|l) \'dU =——— . 
Uid) d{d+l) 



(107) 



By means of our results, it is straightforward to analytically confirm this relation and to 
find the general solution for J \ (1| f/|l) l^dU for arbitrary p G N. In parameterized 
form we have (llf/cll) = ^^^^'^ Y['i=2'^'^^i'^^,n), and accordingly | (1| f/c |1) |^ = 
nn=2 '^^^^(•^i.™) '^^^y depends on \i^2, ■ ■ ■ , ^i,d- Due to this, the integral simplifies as 



follows (see also Appendix Appendix C) 



I (l|f/|l) |W , 

U{d) 

2tt {k>l) i>tt/2 {k<l) d. 
f7r/2 

7r/2 />tt/2 d 



(108) 
(109) 



r ncosP(Ai,„)2(n- l)sin(Ai,Jcos2(--i)-i(Ai,„)c/Ai,2---f^Ai,rf , 
. . . JJ 2{n - 1) sin(Ai,„) cos2("-i)-i+P(Ai,„)dAi,2 ■ ■ ■ d\, . (110) 

7Ai,2=0 •''Ai,d=0„^2 



This integral can easily be solved, i.e. 

d -,2(n-l)+p 



\{2{n-l] 



cos 



(Ai,„) 



_n=2 
d 



2{n - 1) +p 



7r/2 



Ai.„=0 



2{n - 1) 



2(n-l 



n=2 



+ P 



(111) 

(112) 



For the special case p = 4 the general solution 



/ i(iif/ii)iw=n-^') 

Juid) n=2 2(^-1 



+ P 



(113) 



simplifies to Y^n=2 2(n+ij — d{d+i) agreement with Ref. [53] • Note that in this way 
we have found a simple necessary criterion for testing if a set of matrices constitutes 
a unitary design. Namely, as there are no distinguished matrix elements and since 
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I (1| f/ |1) = (1| f/ (1| U* it holds: A set of unitaries {Ui}N is a unitary t-design 
only if 

N d 

Y.wum\i)?' = \{^^^, (114) 

^-^ n — 1 + t 

i=l n=2 

for all k,l E {1, . . . , d}, where Wi is the weighting of Ui {wi = for unweighted designs). 
Further criteria can be constructed analogously. 

6.2. Example 2 

It is known that bilateral twirling j U ® UpW ® WdU of a state p on ® results 
in a Werner state [20l [2T| [22[ l3l] which has the form 

i + /3(Et=iKi')0'^l) 

= — d(dW) — ' ^'''^ 

where —1 < /3 < 1. Using a symbolic computation software, we explicitly and 
analytically twirled the maximally entangled state 1^) = X]f=i K) ® K) dimensions 
c? = 2, 3, 4, 5 utilizing Theorem [1] and [3] 



/.27r {fc>0 ptt/2 (fc</) 

/ / t/c®f/c|^) (^|t/j8)f/Jrff/rf , (116) 



and alternatively with Theorem |2] and H] 

"271 (fc=/) /-TT (fc>Z) /•7r/2 {k<l) 



Uc(S)Uc\^){'^\Ul,(g)Ul,dUd. (117) 

In all cases, this yielded a Werner state with /3 = 1, demonstrating once more the 
operationality and validity of our results. Note that our results also enable analytical 
twirling of multipartite qudit states. 

6.3. Example 3 

The entanglement of a bipartite qudit state lip) G 'Ha^'Hb = C'^^C'^ can be quantified 

2nj.\\ — d 



via the (normalized) concurrence [551 ESI EZl EH] determined by C (|^)) = -73t(1 — 



d-l " 

Tr(p^)), where is the reduced density matrix pB = Tta{\4') {"ipD- An interesting 
property of a quantum system is the a priori entanglement, i.e. a characteristic such as 
the generic probability that a state is entangled or the average amount of entanglement 
over all states in dependence of the dimension d. Here, we focus on the average 



{c')= [ c'iumdu, 

Ju{d'2) 



(118) 
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where I?/') is an arbitrary state of ® C^. If we express = U {ip) {tpl W as 

W) (V''! = Yli j=i^i^*j 01 obtain the reduced density matrix 

d /d-1 \ 

P^=^ [ ^'^i+ndU*+nd I K) 01 • (119) 
i,j=l \n=0 J 

To solve the integral /^(1-Tr(p|))rff/ we can associate each Ui with a matrix element 
of a c?^ X d"^ unitary matrix. In this way, the average (C^) reduces to a sum of integrals 
over the polynomials I4 = J \ {k\U \l) \'^dU and I2.2 = J \ {k\U \l) p| (m| U \n) \^dU where 
{k\U \ l) and {m\U \n) denote distinct matrix elements. Now, taking account of flllQp . 
it is a simple combinatorial problem to show that in j Tr(p^)(if/ the term 1^ appears (i^ 
times and 12,2 appears 2((i — l)(i^ times. Hence, 

{C') = -^{l - d'h - 2{d - l)d'l2,2) . (120) 
d — 1 

From (11071) we already know that 1^ = ^2(J_|_i) for integrals over U{d'^). It remains to 
determine 12,2 for two arbitrary but distinct matrix elements {k\ U \l) and (m| U \n). By 
computing 



'2,2 



/ \{l\Uc\l)?\{d'\Uc\l)?dU,. , (121) 



analogously to (fTUH]) - (fTT^ one finds with \ {d'^\Uc\l)\^ = sin2(Ai,d2) that 12,2 = 
^2(^2 = 1-^4- Consequently, the average entanglement of a bipartite qudit system is 

This result is graphically depicted in Figure [TJ Physically interpreted it means that the 
higher-dimensional the system is, the more likely it becomes to obtain a highly entangled 
state when picking a pure state of "H = ® C'^ at random. 



7. Summary 

In this paper we adopted the concept of the composite parameterization of the unitary 
group U{d) to the special unitary group SlA{d). We showed that both parameterizations 
can be used equivalently to describe orthonormal vectors and subspaces with the minimal 
number of parameters. The introduced parameterizations are completely factorized and 
therefore beneficial for numerical optimizations. We also determined the infinitesimal 
volume element in terms of the introduced parameters. We derived a general formula 
of the normalized Haar measure for both the unitary U{d) and the special unitary 
group SU{d) of arbitrary dimension. The found expressions give theoretical insights 
into the differential structure of U{d) and SU{d). Moreover, the Haar measure plays 
an important role in all kinds of unbiased randomizations. It was stressed that our 
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d 



Figure 1. The average entanglement (I122p of a bipartite qudit system measured by 
the (squared) concurrence C^(|?A)) for the dimensions d = 2, ... ,12. The average (C^) 
increases with the size of the Hilbert space Ha = (8) C*. 



results also constitute a framework for computing high-order group integrals. By means 
of our approach, we analytically solved several exemplary integrals and found that the 
solutions are in agreement with the literature. As integrals over unitary groups appear 
in various fields from particle physics to quantum optics to quantum information, it is 
to be expected that our results will find several interesting applications. 
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Appendix A. Some explicit expressions for U{d) 
d = 2 : 

Uc = exp(iP2A2,i)exp(iYi,2Ai,2)exp(iPiAi,i)exp(iP2A2,2) 
1 



dU, 



47^3 



sin Ai 2 cos Ai 2 d\i 2d\2 idXi id\2 : 



c 



dU. 



exp (iP2A2,i) exp (iYi,2Ai,2) exp (iPsAg^i) exp (iYi,3Ai,3) exp (iP3A3_2) 
X exp (iF2,3A2,3) exp (iPiAi^i) exp (zP2A2,2) exp (zP3A3,3) , 
1 



4vr6 



sin Ai 2 cos Ai 2 sin Ai 3 cos Ai 3 sin A2 3 cos A 



2,3 



X (iAi,2(iA2,i(iAi,3ciA3,i(iA2,3<iA3,2(iAi,i(iA2,2<iA3,3 
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d = 4: 

Uc = exp {iP2\2,i) exp (iYi,2Ai,2) exp (iPsAs.i) exp (iYi^sXi^^) exp (iPiX^^i) 

X exp (zYi,4Ai,4) exp (iP3A3,2) exp (^^2,3^2,3) exp (zP4A4,2) exp (iF2,4A2,4) 

X exp (iP4A4,3) exp (iY3,4A3,4) exp (iPiAi,i) exp (iP2A2,2) exp (iPaAs^g) exp (iP4A4,4) 
3 

dU^ = — — sin Ai 2 cos Ai 2 sin Ai 3 cos^ Ai 3 sin Ai 4 cos^ Ai 4 sin A2 3 cos A2 3 
4^iu . , . , . , , , 

X sin A2,4 cos^ A2,4 sin A3,4 cos A3_4 (iAi,2'^A2,ic?Ai,3(iA3,i(iAi,4(iA4,i 

X (iA2,3(iA3^2C^A2,4C^A4^2'^A3^4(iA4^3Ai^l(iA2,2'^A3^3(iA4^4 

Appendix B. Some explicit expressions for SU{d) 
d = 2: 

Uc = exp(zZi,2A2,i)exp(zYi,2Ai,2)exp(iZi,2Ai,i) , 
dU2 = sin Ai 2 cos Ai 2 dXi 2d\2 id^i 1 



^2 



d 



Uc = exp(iZi,2A2,i)exp(iFi,2Ai,2)exp(iZi,3A3,i)exp(2Fi_3Ai,3) 

X exp (z^2,3A3,2) exp (zF2,3A2,3) exp (iZi-^Xi^i) exp (z^2,3A2,2) 

4 n 

(if/3 = — sin Ai^2 cos Ai_2 sin Ai^3 cos Ai_3 sin A2,3 cos A2,, 



3 



d 



X (iAi,2(iA2,iciAi,3(iA3,i(iA2,3'iA3,2(iAi,i(iA2,2 



Uc = exp (i^i,2A2,i) exp (iYi,2Ai,2) exp (zZi,3A3,i) exp (zYi,3Ai,3) exp (zZi,4A4,i) 

X exp (iFi,4Ai,4) exp (iZ2,3A3,2) exp (iF2,3A2,3) exp (i^2,4A4,2) exp (iF2,4A2,4) 
X exp (zZ3,4A4,3) exp (zKj,4A3,4) exp (iZi^^Xi^i) exp (zZ2,4A2,2) exp (zZ3,4A3,3) 
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dUi = — sin Ai,2 cos Ai,2 sin Ai,3 cos'^ Ai,3 sin Ai,4 cos^ Ai,4 sin A2,3 cos A 



TT 



2,3 



X sin A2,4 cos^ A2,4 sin A3,4 COS A3,4 (iAi,2'iA2,i(iAi,3(iA3,i(iAi,4(iA4,i 

X liA2,3liA3_2'iA2,4'iA4,2liA3,4(iA4_3Ai_l(iA2,2liA3_3 

Appendix C. Differential matrix representation for U{d) 

Adopting the matrix representation [A^^^] introduced in Section [XT] to differentials, the 
normalized Haar measure may be written as 

d 

dUd = n ^rn,n , (C.l) 
m,n=l 
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with 



r d\i 



[A, 



27r 
2-K 



-c?(cos^(Ai,2)) 



L 2-K 



d\d,d- 

27r 



-rf(cos2('^-i)(AM)) 



-rf(cos2(Ad-i,d)) 

d\d,d 
2-K 



(C.2) 



This notation is useful for the construction of a normahzed Haar measure for problems 
which are independent of certain parameters \k,i (see discussion in Section 13.11 and 
Ref. [I], as well as Example 1). Since in such cases the dependence on A^,; and dX^^i can 
be removed from the Haar measure, one can replace the corresponding entry A^,/ by 
a constant. If one sets ^ = 1 then fIC.ip preserves the normalization of the reduced 
Haar measure. Moreover, this notation could lead to a better understanding of the 
differential geometry of U{d). 



Appendix D. Differential matrix representation for SU{d) 
Analogously, the normalized Haar measure on SU{d) may be written as 



with 



[ Am,,n] 



d\iA 

2-K 
d\2A 



d.>^d,l 
7r 



dU, 



-rf(cOs2(Ai,2)) 



n 



d II ^rn,n ; 

m,n=l 



(D.l) 



2k 

<iXd,d-l 



-rf(cos2('^-i)(AM)) 

-d{cos'^{Xd-i,d)) 
1 



(D.2) 
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